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£ Abstract 

It is generally argued that material classes with inversion symmetry do not produce bulk dipole related 

£Nj second harmonic generation (SHG). So, SHG is then either ascribed to surface effects or bulk related electric 






quadrupole or magnetic dipole effects. Using symmetry and ab-initio potentials we show analytically that 

due to the fact of the decaying harmonic electric field certain diamond crystal orientations, as e.g. Si(ll 1), 

On 

produce a bulk dipole SHG response. For fee and bec lattices with a single atom basis, i.e. for the most 

important metals, however, SHG can purely arise due to the disturbance induced by the surface. Finally we 

propose an experiment, exploiting the different dispersion for the fundamental as well as frequency doubled 

O , radiation to determine this effect. 
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8 The recent high interest in surface science techniques capable to probe surfaces in-situ has re- 

9 newed the interest in linear and nonlinear optical techniques such as second harmonic generation 

10 (SHG). This interest has been accompanied by theoretical efforts to model the measured data, 

11 interpreting them either in a classical picture or taking into account quantum mechanics by in- 

12 volving transitions between initial and final states.^ Especially for SHG the interpretation and 

13 the origin of the data is controversial. Some authors^ claim that part of the SHG response arises 

14 from the surface and bulk quadrupoles or from magnetic dipole effects, whereas others^S are 

15 mainly considering surface contributions. To our knowledge, however, all the literature has up to 
is now neglected bulk dipole transitions in materials with inversion symmetry^^ using the follow- 
17 ing argument: Materials with inversion symmetry are described by potentials with even powers of 
is the coordinates (i.e. r 2 , r 4 , etc.). Because a second harmonic contribution providing response to 

19 an incident field oscillating with frequency u could origin, due to inversion symmetry, only from 

20 forces which are quadratic, or potentials which are cubic in a component of r, this effect has been 

21 directly excluded. Analytically, the quantum mechanical expression for each tensor element of the 

22 second order susceptibility ii^ 
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23 In equ. ([TJ) i, j, fc denote x, y, z, N the density of oscillators, and each sum runs over all inter- 

24 mediate states n, n' . Therefore for fee and bec lattices, as e.g. for Cu or Fe with a single atom 

25 basis, a symmetric potential and therefore well defined wavefunction parity, equ. ([TJ) yields zero, 

26 because then product of uneven number of excitations cannot close the loop back to the original 

27 state (for details see equ. 8.40 in reference 15 ). 

28 To our knowledge it has never been discussed in the context of SHG that the potential of atoms 

29 in the diamond structure is, due to their tetrahedral bonding, intrinsically inversion asymmetric. 

30 This asymmetry can be observed by using a pictorial representation for the tetrahedral hybridized 

31 sp 3 wavefunctions (Figjl^), respectively by displaying the ab initio potential along the <1 1 1> di- 

32 rection (Fig[TJ)), or by plotting the potential within in the (111) plane as equi-contour plot (FigjTJ;). 

33 The potential has been calculated with the density functional theory implementation of VASP 5.2 

34 using generalized gradient approximation (GGA). The interaction between the four electrons and 

35 the Si ion is described by projector augmented wave techniquePl The constructed potential de- 



36 scribing a fourfold ionized Si and four electrons is rather strong (~ — 120eV). The calculated 

37 electron energies, however, are of the order of a few eV, and therefore the electrons "feel" the 

38 asymmetry of the potential. Including the exchange correlation part makes the potential landscape 

39 even more asymmetric. 
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FIG. 1: a) Plot of the probability density for the sp3 hybridized wavefunctions for Si(lll), b) ab-inito 
potential for Si (111) along the [111] direction, c) Contour plot of the ab-inito potential for Si in the (111) 
plane. The red dots represent the positions of the Si atoms. 

40 Plotting the ab-initio potential through another plane, e.g. parallel to a Si(100) surface yields 

41 an inversion symmetric potential and so an electric field vector within this plane cannot excite 

42 SHG. However, because the polarization of the electric field can be in an arbitrary direction to the 

43 bonds, each Si atom will produce for an arbitrary polarization SHG response. Also for the case 

44 of Si(lll), usually the SHG response will be minor, because two neighboring Si atoms form a 

45 symmetric entity, if there is no perturbation. So, in all static or long wavelength limit effects the 

46 SHG response of two neighboring Si atoms cancels, as well as for certain symmetries between the 

47 perturbing electric field and the crystallographic arrangement (e.g. normal incidence on Si(l 10)or 

48 Si(100). 

49 In the general case, however, when the amplitude of the electric field is different on both atom 

50 sites, the two responses will not add up fully destructively and produce in the far field a mea- 

51 surable signal. A strength variation of the electric field could either occur through strong (asym- 

52 metric) focusing, or which is discussed beneath, by the intrinsic absorption of the fundamental 

53 wave, when it propagates into the material. The effect of a varying amplitude of the macroscopic 

54 electric field can be either described by employing spatially dispersive (non-local) models, which 

55 means that the dielectric function depends on the wave-vector and the frequency as in refPj, or, 



56 as done recently, to derive within a classical bond model that spatial dispersion is proportional 

57 to the gradient of the macroscopic electric field along the bond 9 . Using bonds instead of the lat- 

58 tice positions disguises the symmetry properties of the underlying lattice and requires the use of 

59 classical mechanics. Therefore we will proceed by explicitly specifying the spatial dependence of 
so the macroscopic electric fields. If a wave impinges perpendicular onto e.g. silicon with a (111) 
61 orientation (see FIG. |2p), the perturbation through the electric field E can be significant. 




FIG. 2: a) The incoming electric field decays along the propagation direction, b) Layers are shown where 
one Si atom (Si(l)) is covalently bound to one from above and to three Si atoms underneath, where one is 
called Si(2). The center of inversion symmetry (CIS) is located in the middle between Si(l) and Si(2). 



62 Si(l 1 1), respectively its atoms, can be described in a simple picture by four equivalent covalent 

63 bonds with a sequence of layers, which are alternatingly bound by one up bond or by three down 

64 bonds to the next layer (see FIG. [2})). For our model it is insignificant, if we imagine the bonds ei- 

65 ther as rigid, such that the electrons can only move along the bond direction 11 or as sp 3 -hybridized 
ee orbitals^H or as extended coherent wavefunctions, describing the probability density for finding 
67 the electron in a given spatial interval. The electronic probability density is the same for positions 
es around Si(l) respectively Si(2), which are at slightly different z- positions. Therefore, the atoms 

69 experience a different field strength, as indicated in FIG. |2^). So, in the antenna picture for the 

70 bulk (111) case the two waves, originating at Si(l) respectively Si(2) do not fully destructively 

71 interfere. The non-linear polarization for a single atom/layer (SL), either Si(l)or Si(2), is derived 



by 



p SL,2. = | J2 xg(2w; w, u)EMEk{") (2) 

73 where xfk stands for the specific coefficient of the nonlinear hyperpolarisibility of a single Si 

74 atom respectively layer. The nonlinear hyperpolarisibility yields for a single Si(lll) tetrahedron 

75 a finite contribution as it has been found before for single sp 3 - orbital 18 . For the case of the later 

76 discussed simplified bond-hyperpolarisibility (SBHM) model 11 Xik can be easily found by setting 

77 in the cited reference the up- and the down hyperbolarisibilities equal. When the fundamental 

78 field penetrates through the bulk, it decays with a complex wavevector k = coh/c, respectively an 

79 absorption coefficient a w /2, given by ct w = ^p 1 , where A is the vacuum wavelength, and n., is the 
so imaginary part of the refractive index n, connected to the dielectric function e(cu) by e(co) = h(ui) 2 . 
si Because the exciting field decays along z due to absorption, the response of the bonds of Si(l) 

82 does not fully cancel the contribution of the lower atom Si(2). It has to be mentioned that even for 

83 the case that the photon energy is smaller than half the bandgap, the occurrence of SHG implies 

84 absorption of the linear wave. For a bond length in Si of ~ 2. 35 A the period d of the slabs is 

85 ~ 3.16v4 and the normal distance A between the (1) and (2) atoms is ~ 0.8A We calculate the 
as difference of the radiated SHG field for one slab just beneath the surface (z = + ) as the coherent 
87 superposition of the Si(l) and Si(2) atom: 

C ( i) = (Re[e ik " z }) 2 e- a " z Re[e ik2 " z }e- a2 " z / 2 

C (2) = (Re[e^( 2 + A )]) 2 e - Q "( 2+A )Re[e^^ +A )]e- a2 -( 2+A )/ 2 (3) 

p***, = fE ^L (2u .^ u)(C(i) _ Ci2))E iu)E M 

jk 

ss where the first equation yields the response of Si(l) at z, the second equation the response 

89 of Si(2) at depth z + A. We explicitly use real fields in the nonlinear process. The exponential 

90 decaying terms in the field have to be squared in the nonlinear process, yielding the absorption 

91 coefficient. Also the backward propagating wave with twice the fundamental frequency is damped, 

92 which is then detected in vacuum or air. The third equation makes explicit use of spatial dispersion, 

93 i.e. the field dependence, and E® denotes the transmitted space independent fundamental electric 

94 field amplitude at (z = + ), having put the decaying terms into (C(i) — C(2)). We note that 

95 positioning the first Si atom at z, the second one at a depth of z + A is implicitly defining the 

96 disturbed surface region by setting the first unperturbed atom as the one with three down bonds. 



97 Extending the calculation to the full bulk is straightforward by splitting the crystal in slabs with 

98 Si atoms at position (1) or (2), numerating only each double layer by n, replacing z — >■ n ■ d and 

99 summing the fields coherently over all slabs starting with n = 1 as the first subsurface slab. The 

100 procedure is similar to the Ewald Oseen Ansatz^for linear optics with the difficulty of using real 

101 fields. The outside detected field / intensity is the coherent superposition of the field originating 

1 i 1 2 

102 from the surface and the field originating from Si(l 1 1) bulk: I = ^ceo\E sur f + Eb u ik\ ■ 

103 Here we do not treat the field produced by the surface E sur f, where specific models on recon- 

104 struction and electronic states would have to be considered. The bulk contribution, however, can 

105 be calculated analytically by using computer algebra systems (producing long terms), and also 
loo numerically. To determine its relative magnitude, we proceed - in order to show the physics more 

107 clearly - by neglecting the phase of the field (~ sign in equ. ([4])). This analytic procedure is approx- 

108 imately correct for the case of rather high absorption coefficients, when the field is over-critic ally 

109 damped: 

oo 

p&tuu. = fEx sL {2u .^ Uj) E (C7 (1|B) - C7 (2|B) )£J( W )E2( W ) * 

jk n=l ,£. 

jk 

no A modified absorption coefficient asHG = a 2w /2 + a w is governing the contribution of the 

in bulk to SHG. For very large, experimentally not achievable, asHG (penetration depth of just a few 

112 As) the formula above yields a vanishing bulk SHG contribution, leaving only the surface effect. 

113 As discussed before, in the general case the surface and the bulk field should be added coherently. 

114 However, already for harmonic excitation close to the Ex transition in Si, asHG is of the order of 

115 ~ 0.02A -1 , equivalent to penetration depths of just a ~ 120 A for the exciting wave and ~ 60 A for 
us the frequency doubled. For asHG = 0.02A -1 the bulk contribution can give rise to a quarter of the 
117 contribution of a single bond for the probably realistic case of setting the first unperturbed atom 
us as the one with three down bonds. If the first unperturbed atom were the one with three upbonds 

119 for the Si(l 1 1) case, the bulk contribution would give rise to three quarters of the contribution of 

120 a single bond. By this calculation it is clear that the bulk contribution cannot be neglected when 

121 using SHG as surface analytical tool and can have the same magnitude as the one arising from the 

122 surface. 

123 For very small absorption coefficients asnc, i- e - almost transparent materials, we find an ad- 

124 ditional condition, similar to Manley-Rowe relations, which provide the conservation of energy: 



125 Using equ. (4) we use an expansion of the denominator and obtain as leading term x d , equal 

126 to x 5Z/ /4- This result cannot be correct, because vanishing adsorption would then imply a rather 

127 strong bulk contribution. The wrong result originates from two sources: a) for very small ab- 

128 sorption coefficients of the harmonic wave, the hyperpolarisibility depends, due to conservation 

129 of energy, on the absorption coefficient a u . Using an ad-hoc Ansatz for the hyperpolarisibility 

130 to be quadratic -with a small coefficient- in a w , produces the required result that for vanishing 

131 absorption of the harmonic the frequency doubled SHG wave vanishes too. (It also vanishes for 

132 linear dependence), b) The second reason comes from the assumption of neglecting the phases. 

133 This assumption is not justified any more; for propagation lengths of the order of 10/mi phase 

134 (mis-)matching (wavevectors k u , fc 2aJ ) has to be taken into account in equ. Q. 

135 In order to compare the effects to experiments, we now employ the simplified bond- 

136 hyperpolarisibility (SBHM) modeP' to calculate the single atom hyperpolarisibility xfjk- m 

137 SBHM the harmonic polarization is given by P 2w = J2 f3 2 jbjbjbj • mEE — % 2 • *EE, with 

3 

138 the bullets denoting outer tensorial products summing over all bonds j. The hyperpolarisibility 

139 (3 2 j denotes the nonlinear hyperpolarisibility of Si bond, pointing towards the tetrahedra corners. 

140 One of the few wavelength dependent data ((p, P) polarization) for a non miscut Si(l 1 1) sample 

141 have been measured by Kravetsky et alP^ at a polar angle of incidence (AOI) of 45°. 

142 Before discussing their experimental data we note that for small AOIs - or normal incidence- 

143 this model gives rather a 6-fold than a 3-fold symmetry. Assuming a bulk contribution this can be 

144 easily explained, because the projection of the field on the 3 bonds is for low AOIs almost identical. 

145 In bulk Si, due to the high refractive index, the harmonic wave propagates almost perpendicular to 

146 the Si(l 1 1) planes and the polarization vector is almost in-plane, also for an external AOI of 45°. 

147 The experimental results given in ref. 20 are replotted in the left column of FIG. [3} 

us Kravetsky et al. measured the same sample (native oxide on Si(lll)) at three different fre- 

149 quencies, without modifying the surface. Experimentally, it turns out that the SHG signal is, as a 

150 function of the azimuthal angle 0° — 360° either threefold, when the absorption length is small (for 

151 a vacuum wavelength of Ao,2w = 266 nm, FIG.[3](a)) or sixfold when the absorption length is huge 

152 (for a vacuum wavelength of Ao,2w = 532 nm, FIG. [3tc)). The fundamental wave has an AOI of 

153 45° in air, within the oxide ~ 29° and within Si the AOI is, due to its large refractive index, ~ 9°, 

154 depending on the wavelength. 

155 We fit all data of ref 20 with essentially one real parameter: This parameter adjusts the relative 

156 contribution of the surface with AOI = 29°, yielding a mainly 3-fold symmetry, and of the 

7 



157 bulk with AOI ~ 9°, giving rise to a 6-fold symmetry. A second real parameter mimicking 

158 the hyperpolarisibility strength (3 2 , is needed to adjust the arbitrary intensity SHG data. In the 

159 right panel of fig. [3] the fits are shown. The SBHM model yields for the up bond a base line 

160 shift (plotted in blue) and the contribution of the three down bonds are plotted in red, green and 

161 brown. The total intensity is plotted in black. For \o i2uj = 266 nm only an AOI = 29° was used, 

162 for Ao,2w = 532 nm a pure bulk contribution was considered and for the intermediate wavelength 

16 3 \o,2u; = 385 nm approximately 3/4 of the (coherently added) signal arise from the surface (AOI 

164 ~ 29°) and 1/4 from bulk. The assignment "surface dipoles" and "bulk dipoles" in the figure stems 

165 purely from the different AOIs. The fit for Ao,2u> — 532 nm shows a difference to the experimental 

166 data. We assign this difference to the neglection of the wavelength dependence of E sur f. We note 

167 that experimental data of Si(lll) as a function of the azimuthal angle, taken in another groupP^ 

168 are almost perfectly matching our predictions. Despite we use only one (bulk-related) parameter, 

169 it is clear to us that the surface contribution will be important for certain surface terminations- and 

170 furthermore frequency dependent 22 . 

171 Without knowing the absolute SHG intensity it is difficult to determine absolute values for the 

172 bulk hyperpolarisibility. Despite there have been proposals in the literature (refsJ 2 ^ 24 *) e.g. using 

173 steps to separate the "pure" bulk and surface contributions, we propose an unambiguous mea- 

174 surement, separating the harmonic and frequency doubled signals in space. The proposal works 

175 for (slightly) dispersive materials and makes use of the different propagation directions, which 

176 occur in the material for the harmonic and frequency doubled (back-) radiation. (This back radi- 

177 ation is clear in the Ewald-Oseen picture of coherent radiators). Just using Huygens principle for 

178 constructing the planes of equal phase for E u and E 2uJ (see FIG. |4j) one could determine exper- 

179 imentally the bulk SHG contribution by cutting a Si crystal (e.g. with a focused ion beam) in a 

180 shape as shown. 

181 Such a geometry allows to let the E 2uJ field pass through the vicinal interface at a different angle 

182 as the harmonic field, separating these two. Furthermore, through the vicinal cut surface, only the 

183 bulk contributions will contribute to the signal, called 2uJ buik- This signal is again well separated 

184 from the one, where surface and bulk components are in, called Q 2u]sur f- Knowing the angle of exit 

185 and the respective Fresnel coefficient then allows to apply numerically or analytically equations [3] 

186 and to determine experimentally the hyperpolarisibility. Two experimental challenges have to be 

187 overcome: 1) High absorption coefficients cushg an d small propagation lengths within the crystal 

188 will render the radiating region to point-like, and a dispersive source can result; 2) for very small 
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FIG. 3: left column: Measured azimuthal dependence as a function of wavelength (a, b, c,). Data replotted 
(with permission) from ref. 20 ; right column: simulated SBHM data of one slab using one (bulk) hyperpo- 
larisibility and two angles of incidence. For Ao,2w = 266 nm only an AOI = 29° (S1O2) was used, for 
the case of Ao,2cj= 532 nm an AOI = 11.47° for the harmonic wave was used. For Ao,2o; = 385nm a 
linear combination of AOI = 29° and AOI = 11.0° is used. The total intensity is plotted using black lines 
and the intensity from the contributing bonds are given in dashed blue (up bond), red (down bond 1), green 
(down bond 2) and brown (down bond 3) 



189 absorption coefficients the absence of phase matching between the harmonic and second harmonic 

190 wave will lead to partially destructive interference. 

191 Summarizing, we show that in the diamond structure dipole allowed bulk SHG response can 

192 exist, if two atoms, establishing a symmetric entity, experience different field strengths and thereby 




FIG. 4: Proposed geometry for measuring unambiguously the dipolar bulk contribution of SHG. For a 
harmonic wave of Ao,^ = 1064 nm (n=3.64) with an AOI = 45° and a frequency doubled Ao,2w = 532 nm 
(n=4.14, penetration depth ~ 1.3^m) and for 6 v i c of 9.8° the SHG wave passes the vicinal surface without 
being refracted (02ujbuik = 0°), but the harmonic is refracted 9 U = 5.8° to the right. So, the frequency 
doubled wave, which contains surface and bulk contributions, is approximately 35° separated from the field 
only containing the bulk contribution. 

193 radiated harmonic field does not fully cancel. The dipole allowed bulk SHG originates then due to 

194 the finite absorption length of the electric field in material. The magnitude of this effect is of the 

195 order of a single bond and can contribute the same signal as e.g. a surface atom. We compare with 

196 available wavelength dependent SHG data, whose azimuthal dependence can be well modeled by 

197 just using one parameter. Finally we propose a measurement scheme capable of unambiguously 

198 separating bulk and surface effects. 
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